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Abstract 

My colleagues and I identified distant red galaxies (DRGs) with J — K s >2.3 in 
the southern Great Observatories Origins Deep Surveys (GOODS-S) field. These 
galaxies reside at z~l-3.5, ((z)~2.2) and based on their ACS (0.4-1 fxm), ISAAC 
(1-2.2 /im), and IRAC (3-8 /mi) photometry, they typically have inferred stellar 
masses A4>lO n A4q. Interestingly, more than 50% of these objects have 24 /im 
flux densities >50 ^Jy. Attributing the IR emission to star-formation implies star- 
formation rates (SFRs) of ~100-1000 A4 yr _1 . As a result, galaxies with M.> 
10 11 M.q have specific SFRs equal to or exceeding the global value at z~1.5-3. In 
contrast, galaxies with Al^lO 11 Mq at z~0.3-0.75 have specific SFRs less than the 
global average, and more than an order of magnitude lower than that for massive 
DRGs at z~1.5-3. Thus, the bulk of star formation in massive galaxies is largely 
complete by z~1.5. The red colors and large inferred stellar masses in the DRGs 
suggest that much of the star formation in these galaxies occurred at redshifts z>5~ 
6. Using model star-formation histories that match the DRG colors and stellar 
masses at z~2-3, and measurements of the UV luminosity density at z>5-6, we 
consider what constraints exist on the stellar initial mass function in the progenitors 
of the massive DRGs at z~2-3. 
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^50% of the stellar mass in ga laxies today may have 



1 Introduction 

Although as much as 

formed during the sh ort time between z~3 and 1 (e.g. Dickinson et al. . 20031 : 
Rudnick et all l2003h . it is still unclear where these stars formed. E.g., one 
hypothesis is that galaxies "downsize" , with massive galaxies forming most of 
their stars in their current configuratio n at tiigh-z, and lower mass galaxies 
continuing to form stars at lower-z (e.g., Cowie et al. . 19991 : Bauer et al. . 20051 : 
De Lucia etTafl 1200.4 ITuneau et all l2005h . Alternatively, stars may form pre- 
dominantly in low-mass galaxies at high-;?, which then merge hierarchically 
over time, slowly assembling into large, massive galaxies at more recent tim es 



e.g.. lBaugh et all Il99«l [Kauffmann fc Charlotl Il99l ICimatti et alll2002h . 
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Massive galaxies with stellar populations older than a few megay ears have re- 



cently been identified and studied with Spitz er as high as 2~6-7 ([Evles et al 



200,4 Van et "all 1200,4 Mobasher et all |2005| ). One question from this is, what 
are the descendants of these z~6 massive galaxies? And, perhaps more im- 
portantly (especially in the context of this meeting), what kind of stellar pop- 
ulations are forming in these galaxies (i.e., are their enough massive stars to 
produce the Lyman-continuum photons necessary to reionize the Universe)? 
It appears that there are too few l uminous quasars at z>5-6 to provide enough 



global radiation for reionization ( Fan et al. . 2002f ). So, the burden now rests 



on early stellar generations to provide enough UV photons for reionization. 
It may be that there are enough small galaxies at z>6 with high UV escape 
fractions for this task ( Yan et al. . 20031 ). or may be the old stellar pop ulations 
j Q galaxies supplied these photons (Panagia, et al.U2005h . There 



m massive z r 



are as-of-yet few constraints on the stellar populations in these galaxies. 



In these proceedings, I discuss recent Spitzer observations at (3-24 jum) of 
massive galaxies at 2~ 1.5-3 in the southern Great Observatories Origins Deep 
(GOODS-S) Field. The Spitzer /MIPS 24 /zm observations provide new con- 
straints on the SFRs in massive galaxies at these epochs. In addition, the 
Spitzer/IKAC observations at 3.6-8 /iin constrain the total mass in old stars 
in galaxies at these redshifts. As I argue below, the constraints on the mass 
in stars from ancient stellar populations now provide crude limits on the form 
of the initial mass function (IMF) in galaxies near the epoch of reionization. 



2 Spitzer Observations of the Massive Galaxies in GOODS— S 

GOODS is a multiwavelength survey of two 10'x 15' fields, one in the northern 
Hubble Deep Field, and one in the southern Chanrda Deep Field. The GOODS 
dat asets include (along with other things) HST / ACS and VLT/ISAAC imag- 



ing ([Giavalisco et al.l . l2004al) . and recent Spitzer imaging (M. Dickinson et al., 



in prep.). We also use data from Spitzer /MIPS 24 Mm in the GOODS-S field 
from time allocated to the MIPS GTOs (|Papovich et all 12004) . 

From the GOODS-S datasets, my colleagues and I ide ntified distant red galax- 

1031) used this color criterion 



ies (DRGs) with J — K s >2.3 mag. iFranx et al. 



to identify galaxies at 2~2-3.5 that in principle have a strong Balmer/4000 A 
break between the J- and f^-bands. This color selection should be sensitive 
to any galaxy dominated by a passively evolving stellar population older than 



~250 Myr at these redshifts (|Franx et all I2003D . For the GOODS-S ISAAC 
images, the J — K s >2.3 mag selection is approximately complete to stellar 
mass A^^IO 11 A4 & for passively evolving galaxies. In the GOODS-S, we find 
153 DRGs to K s <23.2 mag, spanning 0.8<^<3.7, with (z)~2 .2. Roughly 15% 



of th ese have X-ray detections in 1 Msec Chandra data ((Alexander et al 
20031 ). suggesting that many of these galaxies harbor strong AGN (see Pa- 
povich et al., 2005; L. Moustakas et al., in prep.). 
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Approximately 50% of the DRGs have Spitzer detections with /„(24/im) >50 /iJy. 
The 24 jum emission at these redshifts probes the mid-IR (~ 5-10 ittm), which 
broad ly correlates with the total IR, L IR =L(8— 1000/im) fe.g.. lCharv fc Elbazl . 



2001). We co nvert the obse r ved 2 4 /zm flux densities to a total IR luminos- 
ity using the Dale fc Helou (20021 IR template spectral energy distributions 
and the measured redshifts. Some scatt er is inherent in this co nverstion, and 
is typically a factor of 2-3 in L IR (e.g., b— L^li g. We therefore 



adopt a 0.5 dex uncertain ty on Ltr, which also i ncludes an uncertainty arising 



For the DRGs at z~1.5-3, 

which if attributed to star- 
i 



from redshift errors (see Papovich et al. . 2005). 
the 24 fim flux densities yield L IR «10 1L5 ~ 13 L Q , 

form ation corresponds to SFRs of ^100-1000 M & yr~ (|Kennicuttl . ll99 
2003j ). Thus, a substantial fraction of massive galaxies at these redshifts are in- 
volved in intense starbursts. A similar conclusion is reached fo r -E?zi\ ~-selected 
galaxies at in the northern GOODS field ( Daddi et al. . 2005| ). 



3 Stellar Populations and SFRs in High— z Massive Galaxies 

We modeled the DRG stellar populations by comparing their ACS, ISAAC, 
and IRAC photometry to a suite of stellar-population synthesis models (Bru- 
zual & Chariot, 2003), varying the age, star-formation history, and dust con- 
tent. We use the model stellar-mass-to-light ratios to estimate the galaxies' 
stellar mass. We first allow for star-formation histories with the SFR parame- 
terized as a decaying exponential with an e-folding time, r, ranging from short 
r's (burst of star-formation) to long r's (constant star-formation). We also 
use models with a two-component star-formation history characterized by a 
passively evolving stellar population formed in a "burst" at Zf orm =oo, summed 
with the exponentially-decaying-SFR model above. The two-component mod- 
els check the effects of discrete bursts. Figure 1 shows examples of fitting these 
models to a DRG with an apparent 4000 A/Balmer break between the J- and 
i^-bands. The two panels in the figure show the best-fitting models for the 
one- and two-component models described above. Although the modeling 
loosely constrains the ages, dust content, and star-formation histories of the 
DRGs, it provides rela tively robust estimates of the galaxies' stellar masses 
f Papovich et al. . 20051 see also Forster-Schreiber et al., 2004). In the DRG 
in figure 1, even though the star-formation histories are quite different, the 
derived stellar mass is nearly identical. Typical uncertainties for the stellar 
masses for the full DRG sample are 0.1-0.3 dex. 



Figure 2 shows the specific SFRs (ty/Ai) derived from the masses and SFRs 
the DRGs, w here the SFRs ar e calculated from the summed UV and Spitzer 
IR emission ( Bell et aD . l2005| ). The figure also sho ws the specific S FRs for 
galaxies at lower redshift in the COMBO-17 survey (jWolf et allEool ). which 
overlaps with the GTO Spitzer imaging at 24 fim. The SFRs for the COMBO- 
17 galaxies are calculated using the MIPS 24 /im imaging and rest-frame UV 
emission in an analogous manner as for the DRGs. Masses for COMBO-17 
galaxies were estimated from their rest-frame M(V) and U — V colors, and 



3 




Fig. 1. Illustration of fitting SED models to one of the GOODS-S DRGs which has 
an indication of both young and old stellar populations. LEFT: the best-fit model 
with the star-formation history parameterized as a monotonic, decaying exponen- 
tial with best-fit e-folding time r=10 9 yr. RIGHT: the best-fit model where the 
star formation history consists of two components, a single burst at z=oo, repre- 
senting previously formed stellar populations, summed with a monotonic, decaying 
exponential with best-fit e-folding time r=10 6 yr, representing any ongoing star 
formation. The dashed curves in the right panel show the contribution from each 
of the two model components. The inset gives the best-fit model parameters. The 
data points show the ACS B^Veoe^^zs^o, ISAAC JHK S , and IRAC 3.6-8.0 /im 
photometry and errors. 

have a typical uncertainty of 0.3 dex (see Bell et al., 2004, 2005). 

Although the sample is biased against galaxies with low stellar masses (and 
thus specific SFRs); the data are sensitive to all galaxies with larger stel- 
lar masses. At z~0.3-0.7, there is an apparent lack of galaxies with high 
specific SFRs and high stellar masses. In contrast, the massive DRGs at 
1.5<£<3 have much higher specific SFRs for galaxies. Quantitatively, the 
DRGs with M>10 u M Q and 1.5<*<3 have ^/A4~0.2-10 Gyr" 1 , with a 
mean value of ~2.4 Gyr -1 (excluding X-ray sources). By z~0.7 (0.4) galaxies 
with .M>10 n M Q have \I//.M~0.1-1 (<0.5) Gyr -1 , an order of magnitude 
lower than for the massive DRGs. We define the integrated specific SFR as 
the ratio of the sum of the SFRs, to the sum of their stellar masses, Aii, 
T=J2i ^i/J2i M-i, summed over all i galaxies. This is essentially just the ratio 
of the SFR density to the stellar mass density for a volume-limited sample of 
galaxies. Figure 3 shows the integrated specific SFRs for DRGs at z~l. 5-3.0 
and COMBO-17 at 2 ~ 0.4 and 0.7 with M > 10 11 M Q . The data point 
for the DRGs includes all objects with A4>10 11 A4q, assuming that 24 jum- 
undetected DRGs have f u (24^m) = 60 /zJy, and excluding objects with X-ray 
detections of IR luminosities of color indicative of AGN (see Papovich et al. 
2005). The upper bound of the error box shows what happens if we include 
those objects with possible AGN. The bottom bound of the error box shows 
what happens if we continue to exclude objects with possible AGN, but assum- 
ing that the 24 /zm-undetected DRGs have no star formation. The integrated 
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Fig. 2. Galaxy stellar mass versus specific SFR. The top panel shows the DRGs (red 
circles, grey triangles denote upper limits, black stars indicate X-ray DRGs) and 
galaxies from the HDF-N (blue squares; see Papovich et al., 2005), restricted to 
1.5<z<3.0. The bottom panels show COMBO-17 galaxy samples at lower-redshift 
(as labeled). Open circles show 24 /um-detected galaxies; small filled circles show 
upper limits. Mean uncertainties are indicated as a function of stellar mass. 

specific SFR in galaxies with Ai>\0 11 Ai Q declines by more than an order of 
magnitude from z~ 1.5-3 to z<0.7. This downward evolution in the specific 
SFRs seems to support the so-called "downsizing" paradigm. Our results in- 
dicate that star-formation in massive galaxies is reduced for z<l as galaxies 
with lower stellar masses have higher specific SFRs. 

Figure 3 also shows the specific SFR integrated over all galaxies (not just the 
most massive); this is the ratio of the cosmic SFR density (from Cole et al., 
2001) to its integral, i.e., Y=p*/ / p*dt. The global integrated specific SFR 
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Fig. 3. Evolution of the integrated specific SFR, i.e., the ratio of the total SFR to the 
total stellar mass (from Papovich et al., 2005). The curves show the expected evolu- 
tion of the ratio of the total SFR to the total galaxy stellar mass densities from an 
empirical fit to the evoluti on of the SFR de nsity (solid lines, t hick line includes cor- 
rection for d ust ext inction; Cole et al. , 2001 ) , and the model of Hernauist &: Springell 



(dashed line: l2003l ). The data points show results for galaxies with M* > 10 11 A4q. 
The filled diamonds show the mean values derived for COMBO-17 galaxies, and 
the filled circle shows the mean value for the DRGs. The box about the DRG point 
shows how the results change based on some various assumptions (see text; Pa- 
povich et al. 2005). The error bars do not include systematic uncertainties in the 
SFRs, which are indicated by the inset error bar. 

declines steadily with decreasing redshift, i.e., there is a decrease in the global 
specific SFR. The evolution in the integrated specific SFR in massive galax- 
ies is accelerated relative to the global value. Galaxies with .M>10 n .M© 
were forming stars at or slightly above the rate integrated over all galaxies 
at z~=1.5-3. In contrast, by z < 1 galaxies with .M>10 n A4 Q have an inte- 
grated specific SFR much lower than the global value. Thus, by z<1.5 massive 
galaxies have formed most of their stellar mass, and lower-mass galaxies dom- 
inate the cosmic SFR density (see further discussion in Papovich et al., 2005). 



4 Old Stellar Populations in Massive Galaxies at 2<z<3: 
Implications for the IMF at Higher Redshift 

Because the DRGs have red colors, the model fits favor stellar-population ages 
of 1-2 Gyr. This implies that the earliest DRG prog enitors formed at z>,5-6 



possibly accounting for much of their stellar mass (jForster-Schreiber et al 
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2004: IPapovich et a l.. 2005). Owing mostly to model degeneracies, the uncer- 



tainties on galaxy ages, dust content, and star-formation histories from the 
stellar-population mo deling can be more than an order of magnitude (e.g., 



Papovich et all 120011 ). Nonetheless, with the model star-formation histories 



and stellar masses of galaxies at z~2-3, we can set som e broad statistical con- 
straints on the SFRs of galaxies at z>5-6. For example. iFerguson et al. (2002) 
analyzed the model star-formation histories of LBGs at z~2-3, and concluded 
that in order for the progenitors of these galaxies to produce enough Lyman- 
continuum photons to reionize the Universe requires that the stellar popu- 
lations form from an IMF heavily weighted tow ard high-mass stars, and/or 
an episodic SFR. iForster-Schreiber et al.l (|2004l ) found that if the DRGs had 
a constant SFR over their lifetimes, then integrated UV-luminosity density 
from the DRG progenitors decreases by a factor of ~2 from z~2 to 6. 



It is interesting to ask what the stellar-population models imply for the lumi- 
nosities of the progenitors of the GOODS-S DRGs at z>3-5. Figure 4 shows 
the past SFR for the DRG in figure 1 for the best-fitting models using the two 
different star-formation histories (see § 3, figure 1). The dashed line shows the 
past UV luminosity for the model with a monotonically evolving SFR. The 
solid line shows the past SFR for the model with two star formation compo- 
nents. One component is a monotonically evolving "burst" , which in the case 
of this figure is a very short, recent burst (causing the spike at z~3). The sec- 
ond component respresents a burst of star formation in the past, averaged over 
the time since the Big Bang. A galaxy with two "bursts" of star-formation 
would simply show two spikes in figure 4. If we consider this model as a proxy 
for the history of all galaxies, then the star-formation history in the figure 
represents the maximal rate of star formation due to stochastic bursts as a 
function of redshift for an individual galaxy (see Ferguson et al. 2002). 



Figure 5 shows the integrated past UV luminosity density (a proxy for the 
SFR) for all the GOODS-S DRGs with 2<z<3 and M > 10 11 M Q . For these 
redshifts and stellar mass the sample is approximate complete (i.e., nearly 
volume limited) for passively evolving galaxies. The figure also s hows intanta- 
neous mea surements of the UV lum inosity density at z~3-6 from lSteidel et aL 



(199sj) and lGiavalisco et "aL ( 2004b| ) using LBGs. Comparing the measured UV 
luminosity density to what we infer from the past output from the DRGs, if 
the massive DRGs at 2<z<3 formed their stellar populations with a mono- 
tonically evolving SFR, then they account for <5% of the z~6 UV luminosity 
density. However, if these galaxies instead form stars in quasi-stochastic bursts 
throughout their lifetime, then they account for ^10-15% of the £~6 UV lumi- 
nosity density. Inverting that statement, roughly 10-15% of the stars formed 
in 2~6 galaxies end up in massive DRGs at 2<z<3. 



We can go one step further and impose the constraint that the stars form- 
ing at 2~6 can not overproduce the stellar mass at z~2-3. In this case the 
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Fig. 4. Past SFR for the DRG at z=2.9 illustrated in figure 1. The dashed line shows 
the past SFR for the model with a monotonically evolving SFR, x I / (t)~exp(— t/r). 
The solid line shows the past SFR of a two-component model. One component is a 
recent "burst", modeled as a decaying exponential with a short age and e-folding 
time. The other component approximates the average star-formation history of a 
"burst" sometime in the past between z=2.9 and z~oo. 

data allow us to place a very loose constraint on the IMF of the old stellar 
populations in the DRGs. If the past IMF in these galaxies is weighted to 
high-mass stars (e.g., with a slope shallower than Salpeter, x<1.35), then one 
gets more UV photons produced for a given amount of formed stellar mass 
(more "bang for the buck"). The rate that the light from a stellar population 
fades goes as ~ t~ a . We approximate the change in the past UV lu- 

minosity for a varying IMF by scaling the past UV luminosity by the change 
in 5-band luminosity. For a Salpeter IMF a = 0.8, whereas a "Flat" IMF 
(with slope x=0.0) has a = 1.4. If the old stellar populations formed with a 
top-hea vy (x=0.0), flat IMF, the n the past UV luminosity density is ~90% of 
the z~6 Gi avalisco et al. (|2004bi) value. If the measured UV luminosity 



den sity is in fact lower, t hen the ste llar mass in the DR Gs can account for 
all ( Bouwens et al. . 2005f ). or more ( Bunker et al. . 2004 ) than what can be 



formed with a x=0.0 IMF. Therefore, if the UV luminosity density is lower 
than the Giavalisco et al., (2004b) value, then the IMF must have a slope 
x>0.0 to avoid the awkward situation that the z>6 galaxies overproduce the 
stellar mass in DRGs with 2<z<3 and _M>10 n M & . 

Such a situation seems extreme. It is possible that the massive DRGs contain 
a large amount of the stars formed at 2~6. Currently, we have only the limited 
constraint that the slope of the IMF of the z~6 star-forming galaxies is x>0.0. 
This may be tightened as the observations of the z>6 UV luminosity density 
improve (or we decrease the uncertainty in the stellar-population modeling). 
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Redshift 

Fig. 5. Integrated UV luminosity density from the past best-fit star-formation his- 
tories of all DRGs with 2<z<3 and .M^IO 11 M@. The dashed line shows the past 
UV luminosity density for the single-component star-formation histories; the solid 
line shows the two-component models. Both of the black lines correspond to stellar 
populations formed with a Salpeter IMF, with slope x=1.35. The red line indi- 
cates the past UV luminosity density for the two-component models, and an IMF 
weighted to high-mass stars, with 2=0.0. Data points are from Steidel et al. (1999; 
diamonds), and Giavalisco et al. (2004; squares). 

Moreover, the constraint will be better if we can measure the UV luminosity 
density at yet higher redshift. For example, Bouwens et al., (2004) suggest 
that the z~7-8 UV luminosity density is lower than at by an order of 
magnitude. In this case, based on figure 5 we may already be able to discard 
an IMF with a very shallow slope. We may yet achieve indirect constraints on 
the IMF of high-redshift galaxies at the epoch of reionization using the old 
stellar populations we find at "low" redshifts, z~2-3. 
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